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ABSTRACT 

We explore the stability of tidal streams to perturbations, motivated by recent claims 
that the clumpy structure of the stellar streams surrounding the globular cluster Palo- 
mar 5 are the result of gravitational instability. We calculate the Jeans length of tidal 
streams by treating them as a thin expanding cylinder of collisionless matter. We also 
find a general relation between the density and the velocity dispersion inside a stream, 
which is used to determine the longitudinal Jeans criterion. Our analytic results are 
checked by following the time evolution of the phase space density within streams 
using numerical simulations. We conclude that tidal streams within our galactic halo 
are stable on all length scales and over all timescales. 

Key words: galaxies: kinematics and dynamics - galaxies: star clusters - methods: 
analytical 



1 INTRODUCTION 

Tidal streams are a widespread phenomenon in astrophysics, 
emerging from star clusters (Grillmair, Freeman, Irwin & 
Quinn 1995), dark matter subhalos (Diemand, Kuhlen, 
Madau, Zemp, Moore, Potter & Stadel 2008) or satellites of 
galaxies (Ibata, Gilmore & Irwin 1994). It is present on all 
scales from galaxy clusters (Calcneo-Roldan, Moore, Bland- 
Hawthorn, Malin & Sadler 2000) down to the the very small- 
est dark matter substructures (Schneider, Krauss & Moore 
2010). The gravitationally unbound material forms spectac- 
ular long streams that trace the past and future orbit of the 
host system. 

A wide class of tidal streams can be treated as colli- 
sionless systems, since they are dominated by stars or dark 
matter particles - the local relaxation time within the stream 
is much longer than the age of the Universe. Some streams 
contain gaseous material and are much more complicated to 
understand. For example, the oldest example of a 'stream' is 
the spectacular Magellanic HI stream, trailing well over 100 
degrees behind the Magellanic Clouds. Initially modelled as 
a tidal mass loss feature from the Large Magellanic Clouds 
(Lin & Lynden-Bell (1977)), an alternative explanation is 
that it resulted from a more complex gravitational inter- 
action between the Large and the Small Magellanic Cloud 
prior to their infall in the Milky Way potential (Besla, Kalli- 
vayalil, Hernquist, van der Marel, Cox & Kere 2010). How- 
ever, it may also be the case that this feature is purely hy- 
drodynamical in origin since it contains no stars (Moore & 
Davis 1994). 

Galaxy mergers often create spectacular tidal tails that 



are somewhat different from the streams we consider in 
this paper. These streams are rapidly and violently created 
and they can contain dwarf galaxies aligned along the tails. 
Barnes & Hernquist (1992) carried out simulations of galaxy 
mergers and found collapsed objects populating the stellar 
tails. Therefore they proposed collisionless collapse as the 
creation mechanism of tidal dwarf galaxies. However, Wet- 
zstein, Naab & Burkert (2007) identified these collapsed ob- 
jects as numerical artefacts due to insufficient resolution. 
They rather found that it's the gaseous part of the streams 
that triggers the collapse which leads to tidal dwarf galaxies. 

The dynamics of tidal streams from star clusters and 
dwarf galaxies in our own halo have been extensively studied 
to constrain the mass and shape of the Galactic potential 
(Johnston, Zhao, Spergel & Hernquist 1999; Law, Majew- 
ski & Johnston 2009), alternative gravity models (Read & 
Moore 2005) as well as the orbital history of the satellites 
(Kallivayalil, van der Marel & Alcock 2006; Lux, Read & 
Lake 2010). However, the detailed evolution of the internal 
phase space structure of streams has received less attention 
(Helmi & White 1999; Eyre & Binney 2010). 

Simulations, as well as observational data, show vari- 
ations in the width and the internal structure of tidal 
streams. Likewise the density along a stream can vary con- 
siderably, the most prominent example being the symmet- 
ric streams originating from the globular cluster Palomar 5 
with its equally spaced density clumps (Odcnkirchen et al. 
2001; Odenkirchen, Grebel, Dehnen, Rix & Cudworth 2002). 
There are different explanations for these clumps such as 
disc shocking (Dehnen, Odenkirchen, Grebel, & Rix 2004), 
effects due to the dark matter substructures (Mayer, Moore, 
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Quinn & Governato 2002; Yoon, Johnston & Hogg 2010) or 
epicyclic motions in the stellar orbits (Kiipper, MacLeod & 
Heggie 2008; Just, Berczik, Petrov & Ernst 2009). 

Another interpretation was given recently by Quillen 
and Comparetta (Quillen & Comparetta 2010; Comparetta 
& Quillen 2010), who argued that clumps in streams are the 
result of longitudinal Jeans instabilities. In their model they 
describe a tidal stream as an extended static cylinder of stars 
and they use the results of Fridman & Polyachenko (1984), 
that infinitely extended cylinders are gravitationally unsta- 
ble. With an estimated relation for the velocity dispersion 
and the linear density in the stream, Quillen & Comparetta 
find a longitudinal Jeans length of several times the stream 
width. Comparing their results to the observations of Palo- 
mar 5, they find agreement between the distance between 
clumps in the streams and their fastest growing mode of the 
gravitational instability. 

However, their model of a static cylinder does not take 
into account the expansion that happens due to the diffiu- 
ence of the stars in the stream. Once in the stream the stars 
are no longer bound to the cluster, their intrinsic disper- 
sion causes the stream to grow along the orbital direction. 
Escaping stars also have an intrinsic dispersion, related to 
the dispersion in the outer cluster region. Another way of 
understanding the expansion is by considering the veloc- 
ity difference between the substructure and the outflowing 
stars, which depends on the tidal radius. Since the tidal ra- 
dius is shrinking with time, stars that leave the cluster at 
later times are slower than stars that left before and this 
leads to the expansion of the stream. In reality the situation 
is even more complicated. The stream length is actually os- 
cillating during one orbit, being stretched at pericenter and 
compressed at apocenter. The linear expansion only acts 
on average over several orbital periods. Therefore for short 
timescales, the periodic oscillating effect must be taken into 
account. 

This paper is structured as follows: In section 2 we con- 
struct a simplified model for a tidal stream and we find a 
relation between the stream density and its velocity disper- 
sion. Sections 3 and 4 are dedicated to the study of the 
stream stability, where we first derive the linearised equa- 
tion of perturbations and then look at the one dimensional 
collapse along the stream direction. In section 5 we take a 
critical look at our model by comparing with the detailed 
dynamics of streams using N-body simulations. The orbital 
oscillation of the stream length and its influence on collapse 
are considered. Finally we give our conclusions in section 6. 

2 MODELLING A TIDAL STREAM 

The general case of a streaming cluster is a problem of many 
particle dynamics that can be solved self-consistently with 
simulations. Analytical statements can be made by consid- 
ering a model with simplifying assumptions. Thereby one 
has to be careful to avoid over-simplification. We model a 
tidal stream as a self gravitating cylinder of collisionless mat- 
ter with an expansion in the direction of the cylinder axis. 
For the cluster as well as for the host we choose isothermal 
spheres so that we can use the simplifying relations 
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Figure 1. Orbits of the star clusters in our simulations. The 
eccentricity is given in terms of the parameter b defined as R = bV 
(with 6 = for a circular orbit and b = 1 for a radial infall). In 
increasing eccentricity: continuous (b = 0), narrow-dotted (b = 
0.14) broad-dotted (b = 0.34), dashed-dotted (b = 0.54) dashed 
(6 = 0.74) and continuous (b = 0.88). 

where m, M and cr ci , a gat are the masses respectively the ve- 
locity dispersions of the cluster and the host. The distances 
rt and R are the tidal radius of the cluster and the orbital 
radius to the host. For an isothermal sphere the correction 
factor c ~ 0.8 (Binney & Tremaine 2008). 

Whilst many systems can be reasonably well described 
by an isothermal potential over the scales of interest, our re- 
sults would not apply to systems orbiting within very differ- 
ent potentials. For example, a star cluster within a constant 
density potential would not even produce streams. However, 
for Palomar 5 and for many of the streams in our Galactic 
halo, or within galaxy clusters, an isothermal potential is a 
good approximation over the range 0.01 — 0.5R v i r iai (Klypin, 
Zhao & Somerville 2002). 

In order to test the basic assumptions of our model we 
perform simulations of a star cluster orbiting with different 
eccentricities within an isothermal host potential. The ve- 
locity dispersion is chosen to be a ct — 4 km/s for the cluster 
and Ogai = 200 km/s for the host. Every simulation starts 
with the star cluster at a radius of 20 kpc and we choose 
different perpendicular initial velocities from 283 km/s for a 
circular orbit to 50 km/s for the most eccentric orbit. The 
global potential is a fixed analytic potential whilst the star 
cluster is modelled using 2 x 10 5 stars, set up in an equilib- 
rium configuration at the starting position (Zemp, Moore, 
Stadel, Carollo & Madau 2008). The evolution is followed 
using the N-body code PKDGRAV (Stadel 2001), adopt- 
ing high precision parameters for the force accuracy. The 
softening length of the star particles is e = 0.005 kpc. All 
simulated orbits are illustrated in Fig. 1. 

There are two mechanism responsible for the stream 
growth, on the one hand, the outflow of matter leaving 
the cluster with a certain velocity difference AV and on 
the other hand the stream expansion due to diffluence of 
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the initial dispersion. The expansion velocity w is given by 
w ~ 2a c i, which corresponds to the difHuence velocity of a 
bunch of particles leaving the cluster at the same time. Using 
conservation of angular momentum L leads to the velocity 
difference AV: 

AV rt 

L = RV sine = (R + r t ){V-AV) smO => =f- ~ ^ ( 2 ) 

V H 

Here we have assumed R^> r t and V 3> AV. In an isother- 
mal potential the value of V must be somewhere between 
Vr = (4/n) 1 ^ 2 a ga i and V c = 2 1 ^ 2 a ga i, the radial and cir- 
cular velocities. Using (1) we therefore obtain AV ~ a c ; as 
well as 

w ~ 2AU. (3) 

Physically this means that all particles belonging to the 
stream at to will be distributed over the entire stream length 
at all time t > to. Or in other words, even if there is no more 
outflow from the cluster, the stream always stays attached 
to the cluster. 

The amount of difHuence can be estimated in the simu- 
lation by marking particles at a certain time to and looking 
at where they are found in the stream at t ^> to. The first 
image of Fig. 2 shows a cluster at apocenter after one orbit 
(195 Myr) with the particles of one stream marked in red. In 
the second image we see the cluster at apocenter after nine 
orbits (1756 Myr) along with the distribution of the parti- 
cles marked before. The particles marked at the early time 
are located throughout the stream at later times, confirming 
our above statement. 

The width of the stream depends on the velocity disper- 
sion a. A particle with an energy excess during the outflow 
will be on an orbit with a slightly different eccentricity and 
will therefore complete a full oscillation within the stream 
during one orbital time T. The radius corresponding to half 
of the stream width is then approximately given by 

rx ~ \oT. (4) 

On the other hand the length of the stream after one orbit 
is simply 



lo ~ wT ~ 2<tT, 



(5) 



assuming the approximate relation a c i ~ a. After one orbit 
a single stream should therefore be about twice as long as 
it is wide. This is the case in all our simulations and can be 
checked in the first image of Fig 2. 

The linear density of a stream is given by the relation 

dm rh m , . 

^ = d7 = 7~2AU' (6) 

Here we have used z ~ w ~ 2AU, what results in an ad- 
ditional factor of two compared to a static stream because 
of the stretching effect of the expansion. The rate of out- 
streaming matter is estimated to be 
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- (Ra — Rp 



G a gal 



R, (7) 



Gagal 

where we have used the relations (1). The outflow of the 
matter is averaged over one orbital period. With the relation 
R = bV , where the parameter b depends on the cluster orbit 
(with b — for a circular orbit and b = 1 for a radial infall), 
the linear density becomes 
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Figure 2. Simulation of an isothermal cluster with an eccentric- 
ity of b = 0.74. The image on the top shows the cluster after one 
orbit where the particles of one stream are marked in red. The 
image on the bottom shows the cluster after nine orbits with the 
distribution of the particles marked before. 
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and the Toomre parameter is then given by 



a 

2Gn 



(8) 



(9) 



The smallest value for the Toomre parameter is therefore 
q ~ 1 which corresponds to a radial orbit. The Toomre pa- 
rameter of relation (9) is four times larger than the one 
obtained by Quillen & Comparetta (2010), the reason be- 
ing a factor of two which comes in at equation (6) as well 
as the averaging of the mass outflow in equation (7). Both 
effects are directly related to the expansion of the stream, 
not considered by Quillen and Comparetta. 

An independent way to calculate the Toomre parame- 
ter is by using the virial theorem for an isothermal sphere, 
truncated at the tidal radius r t : 



2 \W\ 

<Tcl = 

m c i 



4nG 



drrp(r)M(r) = — — 
Gm c 



2 

°~cl 



Gm c 
4r t 



Using the approximation a c i ~ a then leads to 



1 = 



2Gm sl 



1 rrid lo 
8 m st r t ' 



(10) 
(11) 

(12) 
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Figure 3. The Toomre parameter as a function of the eccen- 
tricity parameter b. The black dots are the measurements from 
the different simulations. The solid gray line corresponds to equa- 
tion (9), while the dotted line is the prediction from Quillen & 
Comparetta (2010). 



For the extreme case of a radial orbit m c ; ~ 2m st and we 
obtain q ~ 1. This means that for all orbits q must be larger 
than one, a result that confirms the relation (9) above. 

The Toomre parameter can also be determined in the 
simulations by measuring the velocity dispersion and the 
linear density. However, it turns out that the dispersion is 
very difficult to quantify accurately because over one or- 
bital period it strongly fluctuates at any Lagrangian point 
(for example, around any star) . This is due to the oscillation 
of stream-length and stream width, which happens because 
the particles in the stream are on nearly free orbits around 
the host. The velocity dispersion therefore is affected by the 
number of stars used in its measurement since that changes 
the region of the stream over which the dispersion is calcu- 
lated. 

In Fig. 3 we plot the Toomre parameter, where the dis- 
persion is measured in the middle of the stream, at apocenter 
after one orbital period and assuming an isotropic distribu- 
tion (looking at the ration of tangential to radial velocity 
dispersions after one orbit we can see that this assump- 
tion is approximately valid). We notice that the simulations 
roughly follow the theoretical prediction which is given by 
the solid gray line. 

Already at that stage of our analysis it becomes clear 
that the expansion has a strong stabilising effect because 
it leads to a significant boost of the Toomre parameter. In 
the next section we will see that the stability of a stream 
is additionally enforced, since the expanding environment 
leads to a damping in the the evolution of perturbations. 



3 PERTURBATIONS IN AN EXPANDING 
CYLINDER 

In order to find a criterion for the stability, we are now mod- 
eling a tidal stream as a non-rotating elongated cylinder of 
collisionless matter that is linearly expanding in the direc- 
tion of its long axis. For the expansion we introduce the 
comoving coordinate s = az with a(t) — at and set z = lo, 
where lo is the stream length after one orbital period T. The 
expansion factor then becomes 



The orbital period is a natural time measure since the out- 
streaming from the cluster into the tails is mainly happen- 
ing during the cluster orbit from apo- to pericenter when 
the tidal radius is shrinking. During the other half of the 
orbit the tidal radius is growing again and there is nearly no 
streaming mass loss. 

An analytical treatment of the stability of an expanding 
cylinder is possible either on scales much smaller or much 
larger than the cylindrical radius. In the former case we 
can treat the fluid as homogeneous and we therefore get the 
usual Jeans length 



X h j = 



(14) 



With the relation (9) as well as the linear density p — nr^p 
we then obtain 



(15) 



A5 
r± 

Since the eccentricity parameter b is always larger than one, 
the Jeans length exceeds the radius of the cylinder and we 
can exclude collapse on scales smaller than r±. 

However there is still the possibility of collapse in the 
longitudinal direction of the cylinder on scales larger than 
r±. This is the second analytically treatable case which 
leads to a very different stability criterion. In order to de- 
termine the behaviour of longitudinal perturbations we are 
now going to derive the equations for the evolution of density 
perturbations. This is usually done by integrating and lin- 
earising the collisionless Boltzmann equation (Peebles 1980). 
Since we are looking at a thin cylinder, we can assume a 
phase-space density of the form 



ft ,\ \ a[pb + pi{z, 



t)]f(p), 



r < r± 
r > r± 



(16) 



Here we have introduced a homogeneous background density 
pb as well as a first order perturbation pi. An integration 
of the phase-space density immediately leads to the stream 
density 



1 lit, ,\ (A*6+A*i) 
P= - / dpf(z,p,t) = ^ 2 — 



+ (17) 



where D — pi/pt is the dimensionless overdensity. 

The evolution of the phase-space density is described by 
the one dimensional collisionless Boltzmann equation with 
expanding coordinate 

It is now straightforward to derive the continuity and the 
momentum equation of the stars in the cylinder. They are 
given by 

dt(l + D) + ±d M [(v)(l + D)]=0, (19) 

d t [a (v) (1 + D)} + + D) + d z [(v 2 ){l + £>)] = 0, (20) 

where 



(v) = 



J Pfdp 
a J /dp' 



(v) 



J P 2 fdp 
a 2 J fdp 



(21) 
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By substituting the derivative of the second equation into 
the first we finally find the equation of perturbation: 

d 2 t D+2^d t D = ^d z [(1 + D)d z $] + ^d 2 z [(1 + D){v 2 )] .(22) 

In order to solve this differential equation we still need 
to know the potential of a cylinder. The simplest assumption 
is to take 



^\r,z,t), 

i knZ 



b W (r,ty Ko 



(23) 



$(r, z, t) = $ l 
$ (1) (r,2,t) = 

(Fridman & Polyachenko 1984), where fcrj is the comoving 
wave number in z-direction. The Poisson equation for the 
zero-order term is simply 

\d_( d$ (0) s 
r dr 



dr 



= 4irG 



Po, 
0. 



r < r± 
r > r± 



(24) 



with the solution $(°)(r) = TvGppr 2 + const. In contrast 
to a self gravitating cylinder, a stream is embedded in the 
dominating potential of the host and the zero order term 
looks different. However, a dependence of the potential in 
the z-direction only comes in as a first order effect due to 
the internal structure of the stream. Therefore we obtain the 
following Poisson equation at first order 

O^w + I 9r $(D _ ^ = ^ ' (25) 



where pi may vary along the axis of the cylinder. Two inde- 
pendent solutions of this homogeneous differential equation 
are the modified Bessel equations of first and second kind 
Io[x] and -KoM- The general inner and outer solution arc 
given by 



^>(r) = AI 
^\r) = A'l 



kpr 
a 

'kpr 
. a 



+ BK 



+ B'K 



kpr 
a 

kpr 
a 



k 2 



where the boundary conditions require B = A' 
With the matching conditions <&< (rx) 
dr® ( < } (r±) = d r ® ( > } (r±) we find 



(26) 



(27) 



0. 



&>>{r±) and 



4irGpia 2 Ko [k r±/a] 
k 2 W[k r ± /a] 



W - 



kp 



{IpKi + hKo). (28) 



We now look at the case of large perturbations in a thin 
stream (kpr±/a << 1). In the asymptotic limit we get 



2Gpp 



2a 2 



+ 7 + log 



kpr_ 



D. 



(29) 



_(fc r±) 2 " V 2a 
The first order potential inside the stream is then given by 

' k r± - 



q a 



+ 7 



D, 



(30) 



where we have used p, — isr\p together with relation (9). 
The Euler constant is 7 = 0.577. 

Using (22) and (30) we obtain a closed set of equa- 
tions for the perturbations D that can now be linearised. We 
therefore set D « 1, as well as (v 2 ){z,t) = a 2 (t) + 0(v 2 ) 
which gives 



a a 2 



q.2*(l) 



kW 



D, 



(31) 
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Figure 4. Dynamics of the perturbation D with respect to the 
scale factor a for the comoving factors kor± = 0.9 (solid), 0.7 
(wide- dashed), 0.5 (narrow- dashed) 0.3 (dotted) and 0.1 (dashed- 
dotted). From top left to bottom right: q = 0.5, 0.75, 1, 2 



In a collisionless cylinder the longitudinal velocity dis- 
persion decreases as a — apa^ 1 (sec eq. 21), while the per- 
pendicular velocity dispersion stays constant. Equation (32) 
can therefore be written as 



D + 2-D 

a 



2 ;2 

a P K p 
qua 3 



(33) 



The perturbation, D, is damped if the right hand side of 
equation (33) is negative. Therefore we can define a Jeans 
length 



A j = Tvr± exp 



<2<> 



+ 7 



(34) 



which is very different from the stability criterion in a homo- 
geneous surrounding (14). The geometry of a thin cylinder 
leads to a Jeans length with an exponential form that guar- 
antees stability up to much larger scales. 

Equation (33) can now be simplified using (4) and tak- 
ing a as variable: 



D"(a) + -D'(a) 
a 



(35) 



4(fc r ± ) 2 
qpa 3 



{log( 



kpr± 
2a 



+ 7 + 



There are two remaining free parameters, namely qp and 
kor±, which describe the eccentricity of the orbit and the 
scale of the perturbation compared to the width of the 
stream. In Fig. 4 we plotted the numerical solutions for dif- 
ferent sets of parameters. For a small Toomre parameter the 
perturbations will become nonlinear and we expect gravita- 
tional collapse to occur. However, for larger q the perturba- 
tions either undergo a damped oscillation or they freeze out 
after an unsubstantial phase of growth. Comparing these re- 
sults with the relation (9) leads to the conclusion that the 
Toomre parameter of a tidal stream is always large enough 
to assure stability in all cases of interest. 

In order to see the effect due to the linear expansion, 
we also look at the case of a static cylinder. The evolution of 
perturbations is then given by equation (33) with a = 1 and 
a = and its behaviour is plotted in Fig. 5. Even for a large 
Toomre parameter q, there are always collapsing modes sup- 
posing an infinitely extended cylinder. This is fundamentally 
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Figure 5. Perturbations D in the case of a static cylinder, where 
o is the number of orbits. The different lines represent the factors 
kr± = 0.9 (solid), 0.7 (wide- dashed), 0.5 (narrow- dashed) 0.3 
(dotted) and 0.1 (dashed- dotted) in static coordinates. From top 
left to bottom right: q= 0.5, 0.75, 1, 2 



different in the expanding case, where all modes are damped 
for a high enough q, leading to stability on all scales. 

Until now we analysed the stability of a stream with 
linear perturbation theory. In the next section we take a 
different look at the stream stability by exploring the lon- 
gitudinal collapse of cylindrical slices. This somehow more 
heuristic approach is not restricted to the linear regime and 
gives an independent analysis of the problem. 



4 SHELL COLLAPSE IN A CYLINDER 

In a one dimensional case of an extended cylinder the spheri- 
cal collapse reduces to the longitudinal collapse of thin slices. 
We therefore consider a homogeneous and infinitely long ex- 
panding cylinder with a top hat perturbation at the time 
U. The stream can then be cut into slices, which evolve at 
constant energy. The energy at a certain distance s is given 
by 



Ei = \v1 + *(a) = \ 



Since the mass Mi evolves as 



GMi(s) 



Mi (a) 



we obtain the energy 
1 



Ei 



2G^(l + <5). 



(36) 



(37) 



(38) 



Slices with a positive total energy will never collapse and 
therefore Ei > is our stability condition. Equation (38) 
then leads to 



8(1 + <$)ai 

go 



(39) 



where we have used the definition of the Toomre parameter 
(9). Slices further away are stable while nearby ones will col- 
lapse. The critical distance below which the stream becomes 
unstable is growing with the square root of time. 

In the picture of shell collapse the velocity dispersion 
is completely ignored, since the diffusion of particles into 
other slices makes the problem much more complicated. We 



will however account for the dispersion by an ad-hoc intro- 
duction of the Jeans length Aj, which guaranties the stream 
stability on small scales. With (39) and (34) we can then 
construct the stability criterion 



go ^ 



+ -2(q Q a- 1 + 1 ) 



(40) 



which is fulfilled at the beginning (a = 1) and may be vi- 
olated at some later times (a > a c ). This means that for 
ti — to all instable slices are below A,/ and therefore all the 
stream is stable. Later on however and depending on qo un- 
stable modes may appear just above Aj. 

Since the Jeans length gives a minimum size for the final 
structure, the initial collapse must start at a scale well above 
this. A calculation of the collapse-time t co u shows however 
that t co u dramatically grows with the distance of the slice. 
Slices only a few times further away than the Jeans length 
already have a t cM that largely exceeds one Hubble time, 
at least for qo > 1. This means that even though there are 
unstable modes in an expanding stream, they will never have 
enough time to grow substantially. Collapse only occures for 
very small values of qo well below the limit given by (9). 

This qualitative picture is in agreement with the results 
plotted in Fig. 4, where a phase of damped oscillation is 
followed by a phase of growth, freezing out at a very low 
level still in the linear regime. 



5 TOWARDS A REALISTIC STREAM 

A realistic treatment of a tidal stream orbiting its host 
galaxy can become very complex, which leads us to consider 
the possibility that our model of an expanding cylinder is an 
over-simplification and therefore we are missing some impor- 
tant dynamics. In the following we treat possible deviations 
to our model and discuss their influence on the stability: 

• In general, the host galaxy is not simply isothermal, but 
can have a triaxial shape that varies with time, and it con- 
tains substructures. The orbit of a cluster is then no longer 
within a plane and it may lose its regularity. The analysis 
of stability effects in such a complex situation is best tack- 
led with full numerical simulations. Nevertheless, there is 
no a-prior reason to believe that one of these effects could 
fundamentally alter the stability criterion. 

• A stream approximately traces the orbit of its cluster and 
is therefore more and more curved the longer it gets. This 
does not correspond to the straight cylinder used in the 
model. However the effect of the bending is rather stabil- 
ising the stream against longitudinal Jeans instabilities and 
can therefore confidently be ignored. 

• A much more severe limitation of our model is the as- 
sumption of a cylindrical form. In reality streams are often 
more sheet-like and their thickness strongly varies during the 
orbital period. The closer a stream approaches the centre of 
the host, the thinner it gets. The reason for this behaviour 
is the form of the isothermal host potential which leads to 
orbits that occupy a narrower real-space volume closer to its 
centre. In Fig. 6 the image of a stream on an eccentric or- 
bit is illustrated. The difference in the stream-width is very 
pronounced and the sheet like structure at apocenter is also 
visible. Even though the variation in the thickness has a ma- 
jor influence on the local stream density, it does not affect 
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Figure 6. Density map of a star cluster with a leading and tailing 
stream after 2 Gyr in an isothermal host potential. The orbit lies 
in the (y,z)-plane and has an eccentricity factor of b = 0.74. The 
high eccentricity leads to strong variations in the stream width. 
While the streams are narrower and denseer at pericentre, they 
become flattend at apocentre with the typical umbrella-like form. 



the longitudinal collapse condition, which only depends on 
the linear density. Incorporating the effect of the flattening 
of the stream is somewhat more difficult because it affects 
the potential (30). However, it is again unlikely that the 
sheet-like structure would have an enhancing effect on the 
collapse since it is stretching the stream which reduces its 
density. 

• As the stream orbits between apocenter and pericenter, 
its length is oscillating, a fact that is not included in our 
model assumptions and may affect the stream stability. In 
fact, the stream only expands linearly on average, its length 
oscillates during one orbit, being stretched at pericenter and 
compressed at apocenter. In Fig. 7 the average distance of 
random points in streams on different orbits are illustrated 
and the orbital oscillation as well as the overall linear ex- 
pansion are clearly visible. These oscillations have an effect 
on the longitudinal perturbations. From peri- to apocenter, 
when the stream-length is shrinking, we are no longer in a 
stable regime and we expect growth. However, this growth 
happens on a timescale longer than the orbital period so 
that perturbations do not have time to collapse. This can be 
shown by approximating the shrinking of the stream with a 
linearly decreasing scale factor of the form 

o(t) = d-(d-l)|t, (41) 

where d — l ma x/lo- The stream length r(t) = a(t)lo now 
runs from l ma x to lo in half of an orbital period. We then 
use the equation of perturbation (33) and replace the time 
variable with the scale factor. The result is 
2 
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as well as the initial conditions D(d) = 0.1 and D'(d) = 0. 
We find growing solutions if the right hand side of the above 
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Figure 7. The evolution of the distance between chosen parti- 
cles in the stream for different simulations with b = 0.14 (full), 
b = 0.34 (dashed), b = 0.54 (dashed-dotted) , b = 0.74 (narrow- 
dotted) and b = 0.88 (broad-dotted). Whilst there is linear growth 
averaged over the orbital motion, the length is oscillating with the 
orbit, and the amplitude of the oscillation is larger for higher ec- 
centricity. 



Figure 8. Growing perturbations D for a shrinking scale factor a 
withkorj_ = 0.9 (solid), 0.7 (wide- dashed), 0.5 (narrow- dashed), 
0.3 (dotted), 0.1 (dashed-dotted). The plots should be read from 
right to left. Top: q = 0.5 with d = 2 (left) and d = 10 (right). 
Bottom: q = 1 with d = 2 (left) and d = 10 (right). 



equation is positive, where the actual value determines the 
growth rate. A large value of d (high eccentricity) gives a 
small growth factor for a long interval of integration, whilst 
a small value (low eccentricity) gives a large growth factor 
for a short interval, the reason being the cP-term in the de- 
nominator of (42). Hence, the actual growth of perturbations 
stays negligibly small in all cases even for a q as low as 0.5 
and the overall stability of our streams is therefore ensured. 
In Fig. 8 we plotted the evolution of the perturbations be- 
tween peri- and apocenter for the case of q = 0.5 and q — 1 
and with d = 2 and d — 10. 

• Because of the longitudinal contraction at apocenter and 
the transversal contraction at pericenter the density and the 
total velocity dispersion are oscillating twice as fast as the 
stream length. This can be observed in Fig. 9, where we plot- 
ted the longitudinal and transversal velocity dispersion of a 
stream with high orbital eccentricity. The doubling of the 
frequency comes from the fact that the stream is longitudi- 
nally compressed at apocenter and transversely compressed 
at pericenter. Fig. 9 can be understood qualitatively by as- 
suming that the particles in the stream are on nearly free 
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Figure 9. Evolution of the velocity dispersion for a high eccen- 
tricity orbit (b = 0.75,). The dispersion parallel to the stream is 
plotted at the top, the one perpendicular to the stream at the bot- 
tom. The grey dahsed curves show the time evolution predicted by 
the model. The vertical lines correspond to the apocenter passage 
of the cluster. 



tions to grow substantially and the oscillation of the stream 
length has therefore no influence on the stability. 

Collisionless stellar or dark matter streams should 
therefore evolve smoothly in time, simply stretching fur- 
ther away from the parent system. The structure observed 
in tidal streams, such as Palomar 5 must have an external 
origin, perhaps disk shocking or encounters with molecular 
clouds or dark matter substructures. 

Our stability analysis could in principle also be ex- 
tended to other systems producing streams. However, sys- 
tems with non spherical shapes and net angular momentum 
are extremely difficult to analyse with analytical methods, 
since the alignment of the interacting objects is important. 
Merging disk galaxies for example produce streams with in- 
ternal structures strongly depending on the initial alignment 
of the disks and on their angular momenta. In such cases, 
high resolution numerical simulations are the indispensable 
tool for a consistent stability analysis. 



epicyclic orbits around the host, which means that the host 
potential is dominating and that the stream particles are not 
feeling each other. Slightly displaced orbits are then crossing 
at apo- and again at pericenter which leads to large peaks 
in the velocity dispersion. 

Our model predicts a longitudinal dispersion that decreases 
on average, an effect that is not clearly visible in the plot 
on the top of Fig. 9. Whilst the minima in the longitudinal 
dispersion seem to decrease as predicted, the maxima are 
growing with time. This growth comes from the fact that the 
particle orbits separate more and more to end up at distinct 
free orbits with the same eccentricity but with a shift in the 
azimuthal angle. The particles are then all crossing at the 
same place leading to a sharp peak in the dispersion. The 
orbital oscillation is also visible in the plot of the transversal 
dispersion at the bottom of Fig. 9. On average however the 
transversal dispersion seems to stay constant as predicted by 
the model. A more detailed study of the stream dispersion 
was done by Helmi & White (1999), who found a similar 
evolution of the dispersion over many more orbital periods. 



6 CONCLUSIONS 

We have studied the gravitational stability of tidal streams 
by modelling them as thin linearly expanding cylinders of 
collisionless matter. Such a model leads to a stability crite- 
rion that has an exponential dependence on the one dimen- 
sional Toomre parameter. We derive a perturbation analysis 
and also use energetic arguments, to show that a cylinder 
with the dispersion, the density and the growth rate of a 
tidal stream is stable for all times. 

We used numerical simulations to test our main approx- 
imations and to study the detailed phase space evolution of 
tidal streams. As a final consistency check, we note that 
none of our simulations show any evidence for gravitational 
instability. 

In reality, a stream is only linearly expanding on aver- 
age, its length is oscillating during one orbit. This leads to a 
time interval between apo- and pericenter, where the scale 
factor shrinks again and the stream is in an unstable regime. 
Nevertheless, this time interval is too short for the perturba- 
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